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Abstract 

We demonstrate that exclusive B decays to P-wave charmonium factorize in the 
non-relativistic hmit provided that colour-octet contributions are taken into ac- 
count, and estimate the branching fractions. Although there are very large un- 
certainties, we find reasonable parameter choices, where the main features of the 
data - large corrections to (naive) factorization and suppression of the Xc2 and he 
final states - are reproduced though the suppression of Xc2 is not as strong as seen 
in the data. Our results also provide an example, where an endpoint divergence 
in hard spectator-scattering factorizes and is absorbed into colour-octet operator 
matrix elements. 



1 Introduction 



Exclusive two-body B meson decays B M1M2 factorize in the heavy quark hmit 
when the final state meson M2, that does not pick up the spectator antiquark from 
the B meson, is light [1-3]. The reason for this is that the quark- antiquark progenitor 
of the "emitted" meson escapes the decay region as an object of small transverse size, 
which remains invisible to soft-gluon interactions with the B Mi transition [4]. This 
argument extends to any colour- singlet object M2 that is small compared to the inverse 
strong interaction scale, 1/A. In particular, it has been suggested [2] that exclusive 
decays to charmonia if, such as -B — > HK , should factorize in the heavy-quark and 
non-relativistic limit, when the charmonium radius l/(mcf) ^ 1/A. This expectation 
has been confirmed for decays to the S'-wave charmonia J/ip and rjc by explicit next- 
to-leading order (NLO) calculations [5-7].^ However, when the formalism was applied 
to P-wave charmonium states [8-13] infrared (IR) divergences appeared that seem to 
violate factorization. 

In this paper we revisit this problem and show that factorization is recovered, if one 
includes the charmonium bound-state scales ntcV, rricv'^ (with v a characteristic velocity of 
the charm quark in the bound state, v <C 1) into the theoretical framework. These scales 
are assumed to be intermediate between the heavy quark masses rrif,, and the strong 
interaction scale A. The divergence structure described above bears resemblance with 
inclusive charmonium decay or production in the colour-singlet model. As is well-known 
the IR divergence problem in the P-wave colour-singlet amplitudes is resolved by the 
introduction of colour-octet operators [14]. We shall see below that a similar mechanism 
is at work for B — >■ HK . However, there is an important difference between P-wave 
charmonium production in inclusive B decay [15,16] and exclusive decays. While in the 
former the charmonium decouples from the inclusive B ^ X transition below the heavy 
quark mass scale, or at least is assumed to in previous treatments, this decoupling takes 
place for exclusive decays only below the scale of the binding energy mcv"^, since gluons 
with this energy can reconnect to the BK system. In the (formal) heavy-quark limit 
this effect is perturbatively calculable since mcv"^ ^ A, and we shall provide numerical 
estimates for this contribution, which has been neglected in previous calculations. In the 
real world, rricv'^ ~ A, and hence reliable calculations appear to be hard to come by for 
decays to charmonium, contrary to decays to light M2. 

In Table 1 we collect the current branching fraction measurements for the decays 
in question [17-33], including the S'-wave final states. By comparing Xci to the 5"- 
waves, we conclude that the P-wave suppression is almost absent. Recalling that in 
naive factorization only the Xci state is produced, the pattern of the P-wave results is 
even more striking, and suggests significant decay amplitudes beyond naive factorization. 
The colour-octet contributions that we identify here may well be the dominant decay 
amplitudes, although they turn out to be hard to calculate for real-world charmonium. 
Nevertheless, it is interesting to see whether something can be said from theory that 

^Some of these papers use a light-cone rather than non-relativistic description of the charmonium, 
but this does not affect the conclusion. 
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J/i, 




XcO 


Xcl 


Xc2 


K 


B- 


10.26 ±0.37 


9.8 ± 1.3 


1.88 ±0.30 


5.01 ±0.37 


< 0.30 


< 0.38 




8.63 ±0.35 


8.7± 1.9 


< 12.4 


4.0 ±0.4 


< 0.41 





Table 1: Averages of -B ^ charmonium ± K branching fraction data in units of 10~ , 
taken from the Heavy Flavour Averaging Group [34]. The original experimental results 
can be found in [17-33]. The limit for he assumes Br {he — > r/c?) = 0.5. 



may help to understand the pattern of experimental data. 

The infrared divergence in hard spectator-scattering is actually an endpoint diver- 
gence in a momentum-fraction convolution integral. Such endpoint divergences prohibit 
hard-scattering factorization of power-suppressed effects in non-leptonic B decays to 
light mesons [1], and of B to light meson form factors even at leading power [35,36]. 
Understanding whether and how such endpoint divergences can be factorized remains a 
major challenge to theory. In B decays to P-wave charmonia the endpoint singularity 
also arises at leading order in the A/rrib expansion, but we shall show that it can be 
factorized into the matrix elements of colour-octet operators. This is of some conceptual 
interest, since it is not known in general how to factorize endpoint divergences. 

Two recent papers [37,38] also discuss factorization of B decays to charmonium. 
These papers deal with the leading order in the non-relativistic velocity expansion ap- 
plicable to decays to S'-wave charmonia, but do not address P-waves. The term "fac- 
torization" there refers to the heavy-quark mass scale or rric, and should thus be 
distinguished from its use here. While it is evident that B decays to P-wave charmonia 
do not factorize at the heavy-quark mass scale due to the infrared divergences men- 
tioned above, our concern is to show that perturbative factorization is recovered when 
rricv'^ 3> A as conjectured in [2]. Corrections to naive factorization for B decays to 
charmonium have also been estimated with light-cone QCD sum rules [39-41], but with 
this method the issue of infrared singularities in the QCD factorization result is not 
addressed. 

2 Operator definitions and tree-level results 

2.1 Effective Hamiltonian and kinematics 

The effective weak-interaction Hamiltonian for the b ccs transition is 

Hos = + C2Q2) (1) 

with 

Ql = {Cibi)v-A{SjCj)v-A, Q2 = {Cibj)v-A{SjCi)v-A, (2) 
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and neglecting the small contributions from the penguin operators. 

The following notation is adopted for the kinematics of the two-body decay process 
B [xcj, hc)K: pb is the momentum of the B meson, w'^ = p^/Ms with = 1; 
q of the charmonium, with = q'^/M^ the charmonium 4-velocity [Mh being the 
charmonium mass) with v"^ = 1; defines the relative momentum of the c quark inside 
the charmonium, so that Pc = q/2 + I, pc = q/2 — I with q ■ I = 0; p is the momentum of 
the kaon. Since kaon mass effects can be neglected in the heavy quark limit, the vector 
n't = p^ / Ek {Ek being the kaon energy in the B rest frame) is light-like. The opposite- 
pointing light-like vector is denoted n'^ with n_ ■ n+ = 2. We also define z = Ami/ml, 
equal to Mfj/M"^ up to corrections of order^ ^/rrih and v"^. From energy-momentum 
conservation 



E 



■ p 



K 



Ml 



Ml 



2M 



B 



1 



Mf 



M, 



H 



(3) 



For any vector a (or Lorentz index fi) the components transverse to n_ and n+ are 
denoted by a± (/i±), for those orthogonal to v we use aj (/^t)- Thus 



1^ n+ ■ a n 
a , H — n_ 



(4) 



2.2 SCET/NRQCD operator definitions 

The bottom and charm quark masses are assumed to be heavy, with mc/rrib fixed in 
the heavy-quark limit. Integrating out the heavy quark mass scales rub, rric leads to 
an effective theory, in which the b quark is static as in heavy-quark effective theory, 
the charm quarks are non-relativistic (in their center-of-mass frame) as in NRQCD and 
the light quarks are coUinear (or soft) as in soft-collinear effective theory (SCET). The 
situation is similar to the corresponding one for charmless B decays [42-44] , except that 
the meson that does not absorb the spectator quark is now described by non-relativistic 
rather than collinear fields. 

To build the decay operators in the effective theory below the heavy quark scale, we 
introduce the non-relativistic bilinears 



n_„, I --j 75 



2^/2\ 2 



Xv 



75 



Xv) 



^Note the dual use oi v: in power-counting estimates v denotes the smah relative velocity of the 
heavy quarks in the charmonium in the rest frame of the charmonium. In kinematical relations v is the 
charmonium velocity vector. 
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■)M/3q(8) 



Dt 7t 



^lrpA 



Xv 



Xv^ 



Xv 



(5) 



of non-relativistic quark [ip) and anti-quark (x) spinor fields.^ It is convenient to use a 
covariant generalization of the NRQCD Lagrangian, where the non-relativistic fields are 
four- component spinors satisfying = ^Jv and '^Xv = —Xv In the charmonium rest 
frame v = (1,0) the non-zero spinor components reduce to the familiar non-relativistic 
two-spinors, and the T components of a contravariant index equal the spatial compo- 
nents. 

Then, adopting the SCET notation defined in [44], we construct the colour- singlet 
operators 



1 



{CW,)i4W^t]pi_Wc]{sn+){l+j5)K 



OC^'P^'^), (6) 



2S+lp.jl)^| 



(7) 



associated with the vertex and spectator-scattering amplitudes that have been considered 
in previous works [8-13]. The operators O'^C^^^^Pj^^) are written in such a way that 
their tree- level matrix elements are proportional to the form factor f^^{Mjj) times the 
derivative of the quarkonium wave function at the origin, since the expression in square 
brackets in (6) is the SCET representation of the full QCD form factor [45]. The effective 
vertices generated by the A- and 5-type operators (in light-cone gauge, where Wc = 1) 
are shown in Figure 1. New and central to the present discussion are the colour-octet 
operators 



Af2S+l ^(8) 



J 



2S+1 c(8)-] 



Af2S+l ^(8) 



J 



(2S+1 ^(8) 



J 



(8) 



According to the assumption that charmonium is a Coulomb bound state, there exist two 
sets of low-energy scales, rricV, rricv'^ ^ A related to the non-relativistic expansion and 



Dt It denotes the symmetric, traceless part of the tensor Dt It- 
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"W 





^B(^2S+lpjl)-j 



Figure 1: Vertices corresponding to the A- and -B-type operators. 



v^mfoA, A, related to the colhnear expansion and the strong-interaction scale. In the case 
of charmless decays the matrix elements of colour-octet operators can be non-zero only 
due to power-suppressed soft-gluon interactions, where soft means momentum of order 
A, thus they can be neglected at leading order in the A/mf, expansion. For charmonium, 
however, the decoupling of gluons with small momentum holds only when the momentum 
is much smaller than rricv'^; gluons with momentum mcv"^ contribute to the octet operator 
matrix elements even at leading order in A/rrib- These contributions are subleading in 
V, but so are the P-wave operators due to the extra derivative in 0(^^^^Pj^''), hence the 
gluon-exchange contribution to the S'-wave octet operators is relevant at leading order 
in the velocity expansion to P-wave charmonium production. 

Concentrating on the terms relevant to P-wave production at leading (non- vanishing) 
order in the v and A/rub expansion, the effective weak-interaction Hamiltonian below 
the heavy quark mass scale is therefore given by 



where we introduced the short-distance coefficients C and s = ■ p s = Mbs. 
2.3 Tree-level matching of A-type operators 

The leading-order matching coefficients are found by comparing the h — c{p(^c{pc)s 
tree amplitude (topology as in the left diagram of Figure 1) computed with the effective 
Hamiltonian (1) to the corresponding amplitude from (9). The result is: 




Co^(^pf)) 





1 



(2Ci) 




(2Ci) , 
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Cf ^) = \ (2Ci) , Cti^sf^) = (2Ci) . (10) 

Since the matrix elements of the colour-octet operators and spectator-scattering are 
both suppressed by a factor of this reproduces the well-known result from naive 
factorization that only the Xci state is produced at leading order. 

2.4 Estimate of the branching fraction 

The leading-order decay amplitude is now given by the expression 

The hadronic matrix element factorizes at leading order in the expansion in A/rrih and 
as according to 

{HC'^'Pj)K\0\''^'P^/^)\B) = (i^llc'Jf (eX)#+(l -75)/^. - ^/ dsCf^'\s) 



{CWc)ii+[W^J]pi_Wc]{sn+){l+j5)K 



rrib 

B) (i7(2^+^Pj)|O(2^+^Py^)|0), (12) 



where the two factors reduce to the QCD B ^ K form factor, 



(^1 



^r^(ew-c)^+(i - i,)K - ^ / dsc\f\s) 



rrib 

{^W,)i4W^ilp^W,]{sn+){l+^,)K 



2 



1^) = 2Mb ( 1 - ^ 1 iTiMfj), (13) 



Ml, 



and the derivative of the charmonium wave function at the origin: 

{K\OCPi'^)\0) = n^-e{OCPi'^)), 
{Xco\OCP!,'^m = {OCp!,'^)), 
{Xci\OCPi'^m = n.-e{OCP!>'^)), 

{Xc2\OCPi'^m = n^_n^_e,, {OCPi'^)), (14) 

with 



Here spin symmetry of the leading non-relativistic interactions has been used to write 
all four matrix elements in terms of {0(^Pq^^)), or equivalently R2i{0), where R2i{r) 
denotes the radial Schrodinger wave function of the n = 2, / = 1 P-wave states, and the 
prime denotes a derivative. 
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Squaring the amplitude, integrating over the two-body phase space, where we neglect 
the kaon mass, and summing over the charmonium polarizations with the help of 



PolSum[H C^+^Pj)] 



E \n--eiX)\'- 
,=-1 



Ml 
Ml 



2M| 
3^ 



H = Xco 
H = he, Xci 

H = Xc2 



(16) 



we obtain the branching fraction 



Ml-M\ 



H 



167rM| 



A|'2S'+lp(l) 







X PolSum[i7(2'5+^P 



2M^(l-^)/nM^)(0(^p(^)))]'. (17) 



Due to (10) this is different from zero at leading order only for the xa state. With 
parameters as given in Section 5, varying the renormalization scale /i between 2 and 
4.8 GeV, and between 1.45 and 1.75 GeV, we find 



Br {B XciK) = (0.13 - 1.06) ■ 10"^ 



(18) 



The largest uncertainty arises from the scale-dependence of the "colour-suppressed" Wil- 
son coefficient C2 + Ci/Nc- The naive factorization prediction is at least a factor of four 
smaller than the experimental result, see Table 1. This analysis shows that the domi- 
nant contribution to the decay amplitude comes, for P-wave final state, from (naively) 
non-factorizable dynamics, such as radiative corrections, spectator-scattering, and the 
colour-octet contributions. 



2.5 Overview of next-to- leading leading order terms 

The leading contribution to (9) comes from the colour-singlet P-wave A-type opera- 
tors O^C^^'^^Pj^''). Not including the fields themselves in the power-counting estimate, 
this contribution is of order a^v, where the factor v arises from the derivatives in the 
non-relativistic P-wave operators (5) and the factor q;° is from the tree-level coefficient 
function. Our intention is to show that radiative corrections can be consistently com- 
puted, when rricv'^ ^ A, so we now discuss the terms arising at order a^f : 

1. One- loop corrections to the short-distance coefficients C^{^^~^-^Pj^^) of the P-wave 
colour-singlet operators ("vertex corrections"): 

/ihardvcrtcx — _£ T/ T/* ('25'+l p(l) N I / /■^Af2S+l pW\\\ /-|q\ 

^B^H{2S+ipj)K — —/K ^cbVcs'-^ [ -fj Jlonc-loop [ -Tj j;|tree, [i-^) 
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where (...) = {H{^^+^Py>)K\...\B) and {0'^C^+^Py>))\tree equals the factorized 
matrix element (12). 

Loop corrections to the matrix elements of these operators need not be considered. 
In particular, corrections to the factorization relation (12) are suppressed by powers 
of A/m^, V or a^, where the two powers of the coupling arise from the vanishing 
colour-projection of one-gluon exchange. 

2. The tree- level matrix element of the B-type P-wave colour-singlet operators 
QB^2S+ipW-^ is of order as, since one factor oi gs is in the definition of the operator 
and another is provided by the coupling of the coUinear gluon to the spectator- 
quark line ("spectator-scattering"). The B-type form factor is of the same order as 
jBK |-^gj^ hence this contribution is also of order agV. The hard spectator-scattering 
amplitude is 

^ha^d^sptor^^ = / rfSC^(2^+^P«)(S)k,ee(0^(^^-^^Py^)(3))k.ee. (20) 

3. The colour-octet operators 0"^(^^^^Sf^) are of order v^, but their tree-level matrix 
elements vanish. The matrix element is non-zero due to soft gluon exchange be- 
tween the charm-quark lines and the 6-quark and light-quark lines, including the 
spectator quark. At order asV, one must include the one-loop matrix element with 
exchange of a gluon of momentum rricv'^: 

/isoft vertex _, ^soft spectator 

—^VcbV*s ^aC^'^^SP)\tree{O^C^~^^SP))\one-g\uon- (21) 

a=\\,±J'=0,l 

Loop corrections to the short- distance coefficients C^C^^^-^Sf,^) need not be con- 
sidered, since the tree-level matrix element vanishes. 

Colour- singlet S'-wave operators contribute only at higher orders due to the van- 
ishing colour projection in the one-gluon exchange contribution to their matrix 
elements. 

3 Short-distance contributions 

In this section we compute the short-distance contributions related to items (19), (20). 
3.1 One- loop short-distance coefficients 

The one-loop correction to the short-distance coefficients of the A-type colour-singlet 
P-wave operators (6) are obtained from matching the QCD diagrams in Figure 2 to the 
expression (19). There is no contribution from the diagram (not shown in the figure) 
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■)pv\ /Xv tpvX /Xv i^v\ / Xv Vii \ /Xv 





K V e 

Figure 2: One-loop corrections to the short-distance coefficients C^i^^^^P^^^) (vertex 
correction) . 



with gluon exchange between the h- and the s-quark hne, since the B — ^ XaK amphtude 
at tree- level involves only the form factor f^^, which our operator definition reproduces 
exactly, but not /(f'^. 

The loop correction is expected to be large, since it comes with the large Wilson 
coefficient Ci, while the tree amplitudes are either zero or colour-suppressed. We write 

mb^/z All Nc 
rA(3p{i)-._ 1 asCp Ci 3 

C.^(3p-)^-;^^|/[3p.], (22) 

where the extra term in ^Pi^"*) arises from the cc vertex correction in the second line 
of Figure 2. The loop functions f[^^~^^Pj] are extracted from the h{mi,w) — *• c{pc)c{p^s 
amplitude expanded to ffist order in the relative momentum /. The expansion is done in 
the integrand to extract the hard momentum region, and the integration is performed 
after expansion. The decomposition into the four angular momentum states is done 
according to the operators in (5). We use dimensional regularization with c? = 4 — 2e for 
both ultraviolet and infrared singularities, and the NDR scheme (naive anti-commuting 
75) for the treatment of 75 and the definition of the weak effective Hamiltonian (1). 
Ultraviolet divergences are subtracted according to the MS prescription. An ultraviolet 
divergence is present only in /['^Pi], since for ^Pi there exists a non-zero tree amplitude, 
while infrared divergences appear for the other P-wave operators but not for ^Pi. This 
can be understood from the fact that the infrared divergences are related to the ^Sq and 
^Si colour-octet matrix elements, and that soft gluons change angular momentum by one 
unit, but do not change spin. 
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Defining the auxiliary function 

fio = — 



16 Lis 




-1 



Lio 



2(1-2=) 



(1-^ 

81n2 2 41n^z 16 , 2-z , z 



In (23) 



(1-^)2 (1-^)2 (1-^)2 1-Z 2 

and tlie ratio z = Ami/ml of heavy quark pole masses, we find 

1 4(2-3z + 2z2) 4(8-8z + ^2^^3) 
/[ A] = -//o (ln[l - .] - ^n) 

_ 4(6 + 16^ - 19^2 + 5z^) _ 2z{8 + 13z - 18z^ + 5z^) 

(1 -Z)(2-Z)2 (l-^)2(2-z)2 

,r3Di . 4(2-7^ + 2^2) 8(4-^-3^2^^3) 

^[ = ^^"f' " " 

_ 4(6 + 22;z - 27^2 ^ 7^3) 2(20 - 73z + 58^2 _ ^g^S) 

z{l-z){2-zy ^ (1-^)2(2-^)2 

.rs^n 1 2^ 4(3-5^ + 2^2) 

82 , 4^2 

In 2 — In z, 



2-z (l-2)(2-z) 
.rsol 2r Az{2-z-z^) 4z(8 - 8z + 3^2 - 



(2-^)2 (2-^; 



_8z{3 + 8z -9z'^ + 2z^) 2z\l6 - 5z - Az^ + z^) 

(l_^)(2-^)2 (i_^)2(2_^)2 (24) 

The infrared 1/e pole, which violates factorization at the heavy quark mass scale, is 
exhibited explicitly in //o- One should note the difference with decays to light mesons 
or S-wave charmonia, where the same hard vertex corrections are infrared finite. The 
singularity arises here as a consequence of the expansion to first order in the relative 
momentum I, and signals that - unsurprisingly - the colour-transparency argument does 
not hold at the heavy-quark mass scale. 

Results for the hard vertex correction, corresponding to the first four diagrams in 
figure 2, have been obtained previously in [8-13], where the infrared divergences were 
noted for the first time. These papers use a gluon mass rather than space-time dimension 
as infrared regulator, while the finite part is given only in parametric form. This makes 
difficult to compare these result with ours, except for the infrared divergent part, where 
we agree, and for the xa case, where the finite contribution is given explicitly [12] and 
there is no infrared regulator dependence, and where we also agree. The fifth diagram 
in Figure 2 was not taken into account in previous papers. 
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Figure 3: Hard spectator-scattering - tree diagrams for the coefficient functions 



3.2 Short-distance spectator scattering 

The tree contribution to the short- distance coefficients of the -B-type colour-singlet 
P-wave operators (7) is obtained from matching the QCD amplitude for the process 
b{mbw) c{pc)c{pc)s{yp)g{yp) shown in Figure 3 to the expression (20). In this way, 
we calculate directly the momentum-space coefficient function 

C^(2^+ipf )(2/) = I dse'y'C''C'+'Pi'^){s), (25) 

where y = 1 — y is the fraction of longitudinal momentum n^p = mb(l — z) carried by 
the hard-coUinear gluon. We find 



Nc rriby/z \ I - zy ) 

Cf(^^i^^)(y) = ^-fl + ^-). (26) 

Nc mb\ I - zy) 

For later purposes it will be useful to express these results as 

C\-s^^pf\y) = ^ [c^i'^^'Pj] + ^^1^^) (27) 

with I/- independent coefficients C^['^^^^Pj], B['^^^^Pj]^ which follow by comparison 
with (26). 

The hard spectator-scattering amplitude (20) requires the evaluation of the convolu- 
tion 

/ d5C^(^^+^Py))(5)|,ee (0^(^^+^Py))(.))k.ee, (28) 
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Figure 4: Hard spectator-scattering - matrix element. 



which is done following [44]. The corresponding diagram is shown in Figure 4. First, as 
in the case of the A-type operators, the tree matrix element factorizes according to 

{H{''+'Pj)K\0^{''+'pf\s)\B) = {K\ — {m,)i^[Wli]p^W,]{sn+){l+^,)K\B) 

mi, 2 

X (i7(2^+iPj)|O(2^+ipy))|0). (29) 

The heavy-to-light matrix element is the same as appears in non-leptonic decay to two 
light mesons resulting in ([44], Eqs. (6), (8)) 

{K\ — {^Wy^[W^z]p^W,]{sn+){l+j,)h^\B) = 
nib 2 

' dr du dye'^'Ji\{f;y,u;)fB(j)B+{uj)fK<f)Kiy)- (30) 



Arrib Jo Jo 

Here (ps+i^^) and (pxiv) ^i^^e the leading- twist light-cone distribution amplitudes of the 
B and the K meson, respectively, while the hard-coUinear "jet function" J\\{f]y,u!) is 
given at leading order by 

J\'[T;y,u) = -5[T-y). (31) 



Introducing 



rfc^^^-L, (32) 
to Xb 



the tree-level matrix elements reads 



(if|^(eVr,)^[W^t,^xW^c](sn+)(l+75)/l^|S)tree 

irasCp IbMb /"^ , j^^- 1 



/ dye^y'^Mniy). (33) 
Jo V 



Nc nibXB Jo y 

Inserting the Fourier representation of the coefficient function into (28), we find the 
amplitude 

^ hard spectator _ TTQ^Cj? 25+1 p ^ I /0/'2S'+l p(l) \ I n\ 

12 



X f'dyC^C'^'P!j'^){y)^. (34) 

nibAB Jo y 

This result has been obtained previously [8-12] by direct evaluation of the spectator- 
scattering amplitude. 

The problematic aspect of this expression is the endpoint divergence of the convo- 
lution integral when y — > 1. The kaon light-cone distribution amplitude behaves as 
(pxiy) ^ y for small y. But contrary to the situation for decays to two light mesons or 
an S-wave charmonium and a kaon, the coefficient function C^{'^^~^^Pj^^){y) contains a 
piece proportional to 1/y resulting in a logarithmically divergent integral. We regularize 
this integral by introducing a cutoff that replaces the upper limit by 1 — /i with ^ 1. 
Using 

/'-"..^^fel^^. r,, tefa)+fk-(l) ^oM (35) 

Jo Jo 

the convolution integral in (34) takes the final form 

r'rf,C-(^--^Pf)(,)M^ = ^|5[^-+^P.]0k'(l)ln/i 
Jo y J\c [ 

+ B[^^-^^Pj] ("dy'tMAp^^ + c'^^^s^^P,] tdy^]. (36) 
JO y Jo y ) 

The regulator-dependent In /i term appears to violate factorization even at leading power 
in the heavy quark expansion. We shall show below, however, that this dependence is 
canceled by a corresponding ultraviolet divergence in the colour-octet matrix elements. 

4 Colour-octet matrix elements 

The important new element in our treatment of B decays to P-wave charmonia are the 
colour-octet contributions (9), (21) to the decay amplitude. We shall now show how to 
compute these matrix elements when the charmonium is a Coulomb bound state and 
demonstrate that the infrared singularities in the vertex correction and hard spectator- 
scattering can be absorbed into a renormalization of the colour-octet matrix elements. 
In the Coulomb limit the octet matrix elements can be computed in perturbation theory, 
and we provide results at order as, corresponding to the accuracy of the short-distance 
terms. In general, and more realistically, the colour-octet matrix elements may be intro- 
duced as new non-perturbative parameters, but their scale dependence still cancels the 
factorization scale dependence of the hard-scattering terms. However, factorization in 
the sense of separating the B K transition from the vacuum to charmonium matrix 
element holds only in the Coulomb limit, since otherwise the octet matrix elements of 
the SCET/NRQCD four-quark operators contain strongly interacting gluon exchanges 
between the charmonium and the BK system. 
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Figure 5: Amplitude exhibiting a quarkonium pole after summing Coulomb gluon ex- 
changes in G^p. 



4.1 Reduction formula for quarkonium matrix elements 

We briefly review the formalism for the calculation of quarkonium matrix elements 

{HX\0\Y), (37) 

where O is some operator, H the quarkonium with momentum q = Mhv in some polar- 
ization state, and X,Y denote arbitrary other particles. 

The quarkonium matrix element is identified in the standard way as part of the 
residue factor of the pole of a suitable Green function. We pick a Green function with two 
external charm quark fields. The quarkonium bound state poles appear after summing 
ladder diagrams with colour- singlet Coulomb (potential) gluon exchange as the bound 
state poles of the Coulomb Green function G^^^ as indicated in Figure 5. For external 
momenta pi,p2 with [pi +^2)^ near 4m^ the bound state poles appear in the integration 
region where qi = rricV + ^1, q-2 = rricV + £2 with E = v ■ {qi + q2) — 2mc = v ■ {ii + £2) 
and iiT, ^2T = —hr = —£t small, of order rrica'^ and mcOs, respectively. A specific 
charmonium state is extracted by performing the following steps: (1) Insert the spectral 
representation of the Green function and use that near a bound state pole 

G^^\k,£; E) V (38) 

, E.„ — E — ie 

where '^nimi'^) is the momentum space wave function of the Schrodinger operator in 
the nlm basis. The three-vectors are introduced by writing four-vectors aj orthogonal 
to i; as ay = L'^j^{v)a'^, where = (0, a) and L^^{v) is the Lorentz boost from the 
quarkonium rest frame to the frame where its momentum is Mhv. (2) To separate the 
degenerate states with different angular momentum and spin, first introduce the spherical 
decomposition 

^nlra{i)=Rnl{£)YU^) (39) 
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defining £ = Tlien perform a Fierz transformation in the Dirac indices of the two 
intermediate charm propagators in Figure 5 to obtain the projection on the spin zero 
and spin 1 components. For / = 1, the J = 0,1,2 states follow from the standard 
Clebsch- Gordon relations. 

The final result can be expressed in terms of the on-shell cc matrix element corre- 
sponding to (37), which we write in the from 



{ciqi)c{q2)X\0\Y) =u,{qi)A{EjT)vc{q2). 



(40) 



Defining the matrices 



A[H] 



75 



-4«(A)?7^ 



75 



H = K 
H = Xco 
H = Xci 

H = Xc2 



(41) 



in Dirac-index space (and diagonal in colour space), the desired quarkonium matrix 
element is 



{HX\0\Y) = J2Mh 



1 



dH 



3 ~ 



y/2Nl J (27r)3 V 47r 



— i?2i(^) tr A[i7] 



1+ 



1- f 



■ (42) 



This is valid in a leading-order treatment of the non-relativistic bound state dynam- 
ics. Beyond this approximation, corrections to the wave-function and trace expression 
are required. Eq. (42) can be used to calculate the first non-vanishing contribution 
to a quarkonium matrix element, and this will be sufficient for the colour-octet terms 
considered below. 

The momentum-space radial wave function follows from the Fourier transform of the 
position-space expression Rni{r)Yirn{yt) and is given by 



Rnl{k) 



drA-Kr'^ji{kr)Rni{r). 



For the case n = 2, / = 1, using the spherical Bessel function 

sin kr — kr cos kr 



ji{kr) 



{kry 



(43) 



(44) 



the integral evaluates to 



5 fu\ ( ■\ 10247r7Bfc 



(45) 
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where -R2i(0) denotes the derivative of the position-space wave function at the origin, 
and 

IB = ^ (46) 

is the inverse Bohr radius of the charmonium. 

As an example of using (42) we calculate (xci|C(3pf ))|0) ^j^^ 0{^Pi^^) given in (5). 
The tree amplitude is 

AiE,ir) = ^[fT,7^]l, (47) 



resulting in 



tr ( MXci] ^ A{E, ir) ^—A =^in.. e(A). (48) 



It follows that 



n„ ■ 6(A) ^2N,pM^^, \l— i-z) R',,{0), (49) 



which is consistent with (14), (15). 



4.2 Soft vertex correction 

We proceed to the calculation of the colour-octet matrix elements {xcjK\0^{^Si^^)\B) 
and {hcK\0^{^SQ'^)\B) {a =||,±). Note that at order as the interactions are spin- 
symmetric, so there is no contribution of the ^5*1^^ CSq"^) operator to the he {Xcj) final 
state. Each matrix element receives contributions from vertex diagrams (first four dia- 
grams in Figure 2) and spectator-scattering diagrams (Figure 3 with gluon attached to 
the undisplayed spectator quark line), except that now the gluon virtuality is small, of 
order {jricv'^Y vertex diagrams. 

We begin by writing down the soft gluon coupling to the cc pair. The leading inter- 
actions of dynamical gluons with momentum of order ttlcv'^ to the heavy charm quarks 
are provided by the (P)NRQCD effective Lagrangian. They read Qgip'^ {x){A^{t^Qi) — x ■ 
E{t, 0))ip{x) together with a similar term for the antiquark field. The contribution from 
the coupling cancels in the sum of the attachments to the c and the c line (or can be 
gauged away), leaving the chromoelectric dipole interaction. The dipole interaction pro- 
vides the additional factor of velocity v which renders the octet S'-wave operator matrix 
element of the same order in v as the singlet P-wave operator. 

The part of the amplitude involving the charm quark lines can now be expressed in 
the form 

A,,{EJt) = / ^T,[''+'sf^]tG'^^\f-r,kT + iT;E + vk)5Uik + iJ), (50) 
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0A(2S+15(8)) 



Figure 6: Soft gluon attachment to the charm quark hues. 



see Figure 6. The various items in this equation are as follows: k denotes the outgoing 
soft gluon momentum; r(j[^'^"^^5'j^''] comes from the cc part of the colour-octet operator 



as given by the contents of square brackets in (5); 6U{k + 
soft gluon interaction vertex 



is the momentum-space 



6U{k + 



-igsT 



B 



d_ 



9t^ ■ 



(51) 



where p and "B" are, respectively, the Lorentz and colour index of the soft gluon; 
G'^^\ij, kj + '(^T] E + V ■ k) is the Coulomb Green function that sums an infinite number 
of gluon exchanges "between" the operator and the soft gluon vertex. Here we need the 
colour-octet Green function, since the cc pair is produced at the operator vertex in a 
colour-octet state. The calculation of the vertex diagrams with the full Coulomb Green 
function is quite involved (see [47] for the calculation of colour-octet inclusive quarkonium 
production matrix elements), but turns out to be unnecessary to good approximation. 
Similar calculations involving colour-octet Coulomb Green functions [48,49] find that the 
numerically largest term arises from the no-gluon exchange term, since every colour-octet 
exchange is suppressed by the small colour factor —l/{2Nc). We therefore simplify (50) 
by approximating 



rric 



obtaining 



E + V ■ k + — + ie] 
rric / 



(52) 



rpB 



(53) 



Attaching the gluon to the bottom quark and the strange quark line and making use 
of (42), we arrive at 



2M, 



H 



^/WJ (27r)3 V 47r 



{2-kY F + ie 
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X tr (^A[H] ^ A,,{Eh, h) ^) 

X t'gs - {mWc) Ml - i,)T^T^K\B), (54) 

\PB-k p-kj 

where we used that k is soft to simphfy the bottom and strange quark propagators, and 
denote the binding energy by 

Eh = Mh- 2mc + ie. (55) 

The colour index "A" is contracted with an index hidden in the definition 
while p and "B" are contracted with the corresponding indices in (53). The order of the 
colour matrices T^T^, which is different for the attachment to the bottom and strange 
line plays no role, since the colour part of the trace evaluates to ^S"^^ . The corresponding 

matrix elements of the operators 0^(^^~^^ Sf^) have 7_l^ instead of /i+ in the {K\...\B) 
matrix element. The absence of a transverse vector implies that the matrix element 
vanishes, so at this order 

{HK\OtC'^'sf^)\B)U,,^ = 0. (56) 

The integrations appearing in (54) can be performed exactly. The loop integral is of 
the general form 

/d^k k i 

(2^ k^P.k{v.k + u;r = (i^ ^'^'^ + ^'"^'^ ■ ^^^^ 
Only Ip enters the final result, and can be calculated directly by introducing Feynman 
parameters. In contrast to the hard vertex correction, this loop integral is infrared 
finite. The sensitivity to very soft gluon momenta is cut off at the scale of the binding 
energy Eb since u = Eb — ij/rric. On the other hand, there is a logarithmic ultraviolet 
divergence, which we regulate dimensionally to be consistent with the calculation of the 
short-distance correction. The angular ^-integral is also easily done using^ 

dne'^e^j = -^fgf. (58) 



The remaining ^-integral is of the type 

2{EH-f/m, 



roo _ 

/ dii''R2i( 
Jo 



(59) 



A + B\n 

with ^-independent constants A, B, and is evaluated using the explicit form (45) of the 
momentum-space radial wave function. Including the colour and spin traces we find 



" Air i\c 

1^^3c(8)md\I _ Ois Cf .3 (1), , (1), 



tree 5 



{XcjK\Of{'Sn\B)U,,^ = -^::-^VU'Pr)Fs. {xcjK\0^{'Pr)\B)U, (60) 
*Note il = = 



An Nc 
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where we have re-expressed the product of the derivative of the wave function at the origin 
and the B K form factor in terms of the factorized matrix element (12) to facilitate the 
comparison with the hard vertex amplitude (19). Eq. (60) contains the spin-dependent 
coefficients V,^{^p[^^) = 1 and Vs^i^P^^^) = {l/(v^2), 0, 1} (for J = 0,1,2), and the 
loop function 



Fsv = — (1 - z + lnz) 

mb{l - zY \ 



1 

- + In At + In 2 - 2AF 



1 TT^X 

+2112(2) - -Ixi^ z + 2\nz\n{l - 2) - 2 In 2 - j , 

= ^ 3 (2V-m.E. + 7B)- • ^''^ 

This includes the ultraviolet divergence explicitly, as well as dependence on the bound 
state parameters 75 (46) and Eh (55) in AF. The colour-octet matrix elements are 
complex, since they contain soft-rescattering phases. However, as can be seen from the 
expression for AF or from (59), in the present one- loop approximation a rescattering 
phase exists only for positive binding energy. 

Consistency of the approach requires that the infrared singularities in the coefficient 
functions of the P-wave colour- singlet operators cancel with the 1/e pole in (61). In this 
case we may interpret the singularities and corresponding //-dependence as factorization 
scale dependence that cancels in the unambiguous sum of the two contributions. Includ- 
ing the tree-level matching coefficients for the colour-octet operators (10) and making 
use of (22), the cancellation condition reads 

C^i 25+lp(l)) _ 25+1^(1)) ^ ^ 25-Mp(l)) ^ ^^^^ (62) 



' c 



Inserting results from (22), (23), (61), one checks that the 1/e terms do indeed cancel. 
In particular, the vanishing of Vsy{'^Pi^^), hence the complete absence of a soft vertex 
contribution, for the Xci case is consistent with the absence of an IR divergence in the 
loop coefficient f[^Pi]. 



4.3 Soft spectator-scattering 

Next we calculate the spectator-scattering contribution to the colour-octet matrix ele- 
ments. We first consider the part of the amplitude shown in Figure 3, before attaching 
the gluon to the spectator quark line. In the present tree approximation, the gluon 
momentum k is the difference between the momentum of the antiquark in the kaon, 
and the spectator-antiquark momentum p? in the B meson. All components of these 
momenta involve factors of A, except for n+ ■ p^ = yn^p = 2yEx- Since rricv'^ ^ A by 
assumption, we may drop all small components (except in the denominator of the gluon 
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propagator, which would be exactly zero) and approximate k by 



n_ 



k^' = n+-k^ (63) 

with ■ k mcv"^, since k is soft. This implies that ?/ ~ u^, so soft-spectator scattering 
corresponds to an endpoint configuration, in which very little momentum is transferred 
to the spectator antiquark. Almost all of the kaon's momentum is carried by the quark 
generated at the b ^ s vertex. The gluon virtuality is given by k"^ = — n+ ■ n_ ■ p? = 
-2EKyuJ. 

The starting expression for the matrix element is 

X i- /; dy /. My) ^ tr (a[..] A.-.iE., ..) 1^) 
X (-z)^?,tr (T^T^) (-^) MstT f ^757^^(1 - 75) ^ ^^+757.1 , (64) 



4 V 4 ^ 

where the second trace arises from the gluon coupling to the spectator quark and the 
projections on the leading-twist light-cone distribution amplitudes of the kaon and the 
B meson. The attachment of the soft gluon to the cc lines, Acc{EH,iT), is the same 
for the vertex and spectator contribution. Substituting (63) into (50) shows that the 
gluon index p must be transverse. The corresponding matrix elements of the operators 
O'^C^^'^^Sf^) have fl^ instead of 7_l^ in the second trace. The trace then vanishes (also 
for the projection on 0b„(ci;), since a must be transverse), so at this order 

{HK\OlC'^'sf)\B)U^,,, = 0. (65) 

Further evaluation of (64) is straightforward: perform the traces; convert the u- 
integral into 1/Xb (32); do the £ angular and then the radial integral. To facilitate the 
comparison with (34), we provide the final result for the partial amplitude 



A7Z'^PslZ^)K = ^V^,K:Ctr^^S]?>) (M|0^(-+^5r)|S)|.pect (66) 
rather than the matrix element itself: 



^ soft spectator _ t/ T/* ^J^^f^ / t7/'2S'+l p n | /-/^ /25+1 p(l) n | n\ 



X 



2^Eh \ , 1B 



+ 



(67) 
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Here i?[^'^+^Pj ] is given by the same expression as defined in (26), (27). Notice that 
there is no soft spectator-scattering contribution to the he final state, as B[^Pi'\ = 0, 
which will be important in the numerical analysis. 

We compare the integral over the kaon distribution amplitude to (34). While the 
integrand there was applicable to y not near 1 and exhibited a logarithmic endpoint 
divergence as |/ — *• 1, the present integrand is appropriate only to 1 — ?/ ~ w^, i.e. in 
the endpoint region. There is no divergence here as y ^ 1. However, for y ^ v'^ the 
integrand has the same logarithmic behaviour / dy cf) f^{y) / y"^ as does the hard-spectator 
contribution for y <^ 1. In (35), (36) we regulated the endpoint divergence in hard 
spectator-scattering by cutting off the y integral above 1 — /x. This corresponds to a hard 
factorization scale in the energy of the gluon that connects to the spectator quark. The 
spectator-scattering contribution to the colour-octet matrix element originates precisely 
from the energy region that was left out in (36), thus the correct interpretation of the y- 
integral in (67) is /q^ dy J^_^^ dy. To combine with (36) we must evaluate the regularized 
version of (67) up to terms of order v'^/fi. This allows us to approximate (paiy) ~ 
—y(f)'j^{l), resulting in 




^ 0^(l)(-ln/i + F,,), (68) 

where 

Fss = - ln(l - ^) + 1 + zvr + AF, (69) 

and AF is given in (61). Comparing (34) to (67), together with (36), the previous 
equation demonstrates that the regulator-dependent In/x terms cancel. We may there- 
fore conclude that the endpoint singularity is hard spectator-scattering does not violate 
factorization, since it can be factorized into the colour-octet matrix elements. 



4.4 Further remarks on the endpoint singularity 

A rigorous understanding of endpoint singularities in convolution integrals would enhance 
the predictivity of QCD factorization approaches for exclusive B decays considerably. It 
would also provide meaning to formal factorization "theorems" derived in soft-coUinear 
effective theory, which generically result in ill-defined convolutions, with exceptions in 
many leading-power applications. Despite several attempts [46,50,51] there is currently 
no satisfactory framework for factorizing endpoint divergences and for associating them 
with well-defined operator matrix elements. 

The calculation of the B XcjK decay amplitudes presented in this paper provides 
the first example, where an endpoint singularity in a hard-scattering convolution integral 
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can be factored consistently (at least at the leading order) into precisely defined objects, 
the colour-octet matrix elements. The example does not quite represent what is required 
for other cases, such as the B ^ n form factor, since in B ^ XcjK the endpoint diver- 
gence arises from factorization at the scale rricv'^, not A.^ Nevertheless, it is worthwhile 
to collect some observations on the structure of the endpoint contribution. 

• The endpoint contribution is proportional to 0^(1), the derivative of the distri- 
bution amplitude at the endpoint.^ This is because the endpoint region is of size 
v"^ rather than A/rrih, hence it is justified to describe the quarks in the kaon by 
coUinear quark fields. In general, we do not expect the distribution amphtudes 
to be relevant in the endpoint region, since one of the quarks does not carry a 
coUinear momentum. 

• The factor 1/y that renders the hard-spectator convolution integral divergent, orig- 
inates from the expansion of the charm propagators. It is therefore not possible to 
associate the endpoint contribution with a matrix element involving only the kaon 
state, as in the case of the light-cone distribution amplitude. Rather it refiects a 
large non- factorizing contribution to the entire process from the scale rricv'^. 

• The endpoint contribution contains a large rescattering phase as seen from (69). 
This observation casts doubt on the correctness of a claim made in [54] that the 
power-suppressed weak annihilation contributions to charmless decays are real in 
first approximation. In fact, applying the prescription of endpoint subtraction used 
in [54] to B ^ XcjK would simply set the soft-spectator scattering contribution 
to zero. This implies an uncanceled subtraction scale dependence proportional to 
In/i, which is also present in the result of [54]. More important to the claim, it 
would miss the soft spectator-rescattering phase. 



5 Estimates of branching ratios 

Before discussing the numerical results, we present an expression for the sum of the 
different contributions to the amplitude: 

AB^Hes,.p^)K = ^ V^i^K: {HC'^'Pj)\OC'^'p!j'^)\0) [Ay + As] , (70) 

where separating the vertex corrections from the spectator-scattering corrections we have 



Ml, 



CA(25+lp(l))^CA(25+lp(l) 



^This aspect is similar to the discussion of the B ^ rjc form factor in [52], which also exhibits a 
calculable endpoint logarithm. 

^As in the expression for the B ^ ir form factor in the heavy quark limit of its light-cone QCD sum 
rule representation [53]. 
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4^ ivT 



2S+lp(l) 



7^/ 2S+1 c{l)^ 



^ ^ vra^C^ fKfsMB 2Ci i^B,2S+ip 



y 



y2 



(71) 



and all quantities have been defined previously. The branching fraction, extending the 
leading-order expression (17), is given by 



Br {B HC'^'Pj)K) =rBY iVctKll' ZJ!^ PolSum[i/(2^+ip, 



X 



167rM| 



(72) 



5.1 Parameters 

The numerical result depends on the couplings Gp = 1.16639 ■ 10~^GeV~^, iKfeVcsi 



40.4 ■ 10-3. The QCD scale is = 225 MeV (MS scheme, five quark flavours) and 
next-to-leading logarithmic running of the strong coupling and Wilson coefficients is 
used. At /i = 3GeV: a,(3GeV) = 0.2503, Ci(3GeV) = 1.105, C2(3GeV) = -0.2366. 
The renormalization scale for these quantities is denoted by /x. However, in the strong 
coupling that multiplies the spectator scattering term we use the intermediate scale 
Hh = GeV fi, and in the expression for the inverse Bohr radius (46) we imagine 
choosing the scale of (or rric) such that 7^ = 500 MeV. The values of the quark 
masses m^, rric will be discussed below. 

The meson masses are Mb = 5.279 GeV, M^^^ = 3.415 GeV, M^^, = 3.511 GeV, 
^Xc2 — 3.556 GeV, Mh^ = 3.525 GeV. The derivative of the wave function at the origin 
can be determined from Xcj decays and takes the value |-R2i(0)P = 0.050 GeV^ [16,55]. 
The kaon and S-meson decay constants are fx = 160 MeV and Jb = 210 MeV, respec- 
tively, the moment of the S-meson distribution amplitude is assumed to take the small 
value Xb = 200 MeV that is favoured by the large rates of colour- suppressed charmless 
B decays [44]. The B ^ K form factor is parameterized following [56] in the form 



0.162 



0.173 



'1 



(Mb, =5.41 GeV). 



(73) 



Many of these parameters have significant theoretical errors, but in view of other un- 
certainties discussed below, they are less relevant, except for the parameter Xb, where 
values twice as large are not ruled out theoretically. Finally, we expand the kaon light- 
cone distribution amplitude into Gegenbauer polynomials 



n=l 



(74) 
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and truncate the expansion at order n = 2. The first two Gegenbauer moments are 
of (2 GeV) = 0.06 ± 0.06 and a§{2 GeV) = 0.2 ± 0.15 [57-62], while using "asymptotic" 
distribution amplitudes amounts to setting the Gegenbauer moments to zero. In terms 
of Gegenbauer moments the expressions appearing in (71) read 

0^(1) = (-6) (l + 3af + 6af), 

h = dy — - — = 3 1 + af + 02 , 
Jo y ^ ' 

While 1 1 is well-behaved, 0^(1) and I2 exhibit a very large sensitivity to the higher 
Gegenbauer moments. Comparing the maximal value of I2 to its asymptotic one, we find 
jmax^jas _ -j? ^2 and it is not clear whether the Gegenbauer expansion is converging at all. 
A consequence of this is that the size of the spectator-scattering amplitude is uncertain 
by a factor of several (including the uncertainty in \b) for Xcj, where B['^^^^Pj] is not 
vanishing. 

5.2 Results 

Given the large ambiguities mentioned above, but also the fact that our calculation relies 
on the unrealistic assumption that charmonium is a Coulomb bound state, we do not 
expect reliable quantitative results for the B — > {xcj, hc)K branching fractions. Instead 
we address the questions 

1) Are large corrections to naive factorization expected theoretically? 

2) Why are the Xc2K and hcK final states suppressed relative to XcoK, XciKl 

which are of interest given observations summarized Table 1. Our calculation results in 
exactly the same decay rates for charged and neutral B decay. Thus, branching fractions 
of pairs of related decays differ only by the lifetime ratio tb^/tb^^- In the following we 
consider only B^ decay using tb^ = 1.53 ■ 10~^^ s. 

We begin by discussing the dependence of the branching fractions on the various 
inputs, when we neglect the spectator-scattering term entirely. The scale-dependence, 
adopting the quark masses values mf, = 4.8 GeV, = 1.4 GeV, is shown in the left 
plot of Figure 7. There is a large cancellation between the tree level and one-loop 
contribution to the coefficient function relevant to the xa final state resulting in a very 
small branching fraction. The branching fractions for the other final states are also 
quite small, not exceeding a few times 10~^ with he and Xc2 being even smaller than the 
other two. At this point we can already conclude that corrections to naive factorization 
are order one effects, providing a positive answer to the first question above, as is in 
fact expected for colour-suppressed decay modes. The final-state dependence might be 
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Figure 7: Branching fractions in units of 10 ^, when the spectator-scattering contribution 
is set to zero: renormahzation scale (left) and charm quark mass (right) dependence. 



similar to the data, but this cannot be the complete story, since the Xco and Xci branching 
fractions fall short of the data by about an order of magnitude. 

The scale dependence of the NLO result remains significant, simply because the 
LO term for Xci is canceled, and there is no LO term for the other final states. The 
scale dependence is exactly of the form (a^Ci)^ for Xco, Xci-, and he and approxi- 
mately so for Xc\- This causes an uncertainty of a factor of 2 when /i is varied be- 
tween 2 GeV and 5 GeV, and larger if one allows smaller scales. However, below about 
1.5 GeV the scale-dependence blows up as seen in the Figure. In the following we fix 
/i = 3 GeV. Results for other choices of /i can be obtained approximately by multiply- 
ing with (asCi)^(/i)/(asCi)^(3 GeV). (This remains true, when the spectator-scattering 
terms are added back.) Other significant parameter dependences arise from the quark 
mass values. The dependence on rUc is more important than the one on m;,, so we fix 

= 4.8 GeV in the following. The charm quark mass dependence of the branching 
fractions, still omitting spectator scattering, is shown in the right plot of Figure 7, from 
which it is seen that the size of the Xco? Xci branching fractions versus Xc2; he reverses 
as rric increases. The charm quark mass here is the pole mass, which is a poorly defined 
quantity in perturbation theory, due to large radiative corrections. Typical values are 
rric = 1.4 . . . 1.7 GeV. There may be good reason to choose larger values here, since the 
charm quark pole mass controls the binding energy = — 2mc, which should be 
negative in the approximation of charmonium as a non-relativistic bound state, and is 
negative in reality when measured relative to the DD threshold. From the Figure it 
is evident that with the NLO vertex correction alone it is not possible to explain the 
experimental data, since the Xco, Xa branching fractions are too small, even allowing for 
theoretical uncertainties in the B ^ K form factor or the charmonium wave function. 

Next we imagine that the branching fractions are given by the spectator-scattering 
term alone. The largest parameter dependences now arise from the charm quark mass. 
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1.4 1.5 1.6 1.7 1.£ 




1.4 1.5 1.6 1.7 1.8 



Figure 8: Charm quark mass dependence of the branching fractions in units of 10 ^ from 
spectator scattering only. Left for asymptotic kaon distribution amphtude (af- = = 
0), right for af = 0.06 and a§ = 0.2. 




26 



m,/GeV 


XcoK' 








1.45 


1.13 


1.31 


0.28 


0.29 


1.50 


1.17 


1.52 


0.37 


0.30 


1.55 


1.25 


1.85 


0.51 


0.31 


1.60 


1.44 


2.36 


0.74 


0.34 


1.65 


1.95 


3.28 


1.14 


0.38 


1.70 


3.70 


5.23 


1.95 


0.44 


1.75 


5.19 


10.31 


3.98 


0.53 



Table 2: HK^ branching fractions in units of 10 ^ for various charm quark mass 

inputs. Asymptotic distributions amphtudes, \b = 200 MeV and = 4.8 are fixed. 



the Gegenbauer moments of the kaon hght-cone distribution amphtude, and the B me- 
son distribution amphtude parameter A^. In Figure 8 we show the mc dependence for 
asymptotic distribution amphtudes (left) and for = 0.06, af^ = 0.2 (right). The 
branching fractions grow rapidly with rric, when the spectator amplitude becomes domi- 
nated by the imaginary part from the colour-octet contributions. We also observe a huge 
dependence on the Gegenbauer moments, confirming the expectation that the expansion 
may be invalid. Even af = 0.2 leads to unacceptably large branching fractions. These 
effects are less pronounced when As is larger, since the spectator-scattering branching 
fraction shown in Figure 8 is proportional to 1/A^. Independent of these uncertainties, 
we always find that spectator-scattering is a small effect for he, and larger for XcO: Xci 
than for the other two final states. 

When we now add both contributions together, including the interference term, we 
obtain the result shown in Figure 9 for asymptotic distribution amplitudes. For rric in 
the range from 1.5 GeV to 1.7 GeV, this suggests the interpretation that the Xco, Xa final 
states are dominated by spectator scattering, more precisely by the spectator-scattering 
contribution to the colour-octet matrix element. The smallness of the he branching frac- 
tion is explained by the absence of such a contribution (at leading order) for this final 
state. The Xc2 case is intermediate with a rapidly rising branching fraction in the interest- 
ing charm-quark mass window. Numerical results for some values of rric are provided in 
Table 2. We emphasize that in addition to the charm-quark mass dependence displayed 
explicitly there are further large theoretical uncertainties related to scale-dependence, 
which shifts all branching fractions uniformly as described above, to mfe-dependence, to 
the Gegenbauer moments, and to Xb- There are some parameter degeneracies that allow 
making af^ and Xb simultaneously larger. In view of these uncertainties the main con- 
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elusion of the numerieal analysis is that there are reasonable regions of parameter space 
{rric ^ 1.65 GeV, small A^, and asymptotic kaon distribution amplitude), where the the- 
oretical calculations in our model for the colour-octet matrix elements are in qualitative 
agreement with the experimental data, namely the existence of large contributions be- 
yond naive factorization, and the suppression of the Xc2K and h^K modes. From Table 2 
we conclude that the small he branching fraction is a robust feature of our results, but 
we find it difficult to explain the strong Xc2 suppression seen in the data compiled in 
Table 1, while maintaining the sizeable xa branching fraction. 

6 Conclusion 

We revisited exclusive B decays to P-wave charmonia motivated by previous studies 
[8-13] of these decays in the QCD factorization framework that reported a violation 
of factorization. In contrast, we find that after accounting for colour-octet operators, 
which, contrary to the case of charmless decays, are not suppressed by A/ nib due to 
the existence of the charmonium binding energy scale, QCD factorization is recovered, 
at least at order as- The infrared divergences found in previous calculations can be 
subtracted consistently into the matrix elements of these operators. This includes the 
endpoint divergence that is found in the unsubtracted coefficient function associated 
with spectator-scattering. Our calculations demonstrate that the endpoint contribution, 
now contained in the colour-octet matrix element, can lead to a large rescattering phase. 
These observations may be of conceptual interest, since it is presently still unclear in the 
general case, whether and how endpoint divergences that often appear in convolutions 
in coUinear factorization formulas can be absorbed into well-defined non-perturbative 
objects and what these objects are. We find it plausible that factorization of B decays to 
P-wave charmonium extends to higher orders in the coupling expansion when rricV^ ^ A, 
in view of the argument presented in [2]; nonetheless, it would be of great interest to 
verify the factorization of endpoint divergences beyond the tree-approximation to the 
hard-scattering sub-graph. 

Previous numerical estimates of the branching fractions to P-waves relied on ad 
hoc treatments of the infrared regulator dependence. In the present framework, this 
is unnecessary, but an estimate of the colour-octet operator matrix elements is needed, 
which may even be the largest contribution to the decay amplitude. To this end we 
adopted a description of charmonium as a Coulomb bound state, which corresponds to 
the formal heavy quark limit. In practice, this limit is probably unreliable, and our 
results do indeed exhibit large theoretical uncertainties. Nevertheless, we find that for 
plausible theoretical inputs it is possible to reproduce qualitatively what we consider 
to be the most interesting features of current experimental data: suppression of the 
Xc2K and hcK final states and amplitudes that must be dominated by terms beyond 
naive factorization, though the suppression of Xc2K is not as strong as observed. An 
interesting avenue to pursue in the future might be to consider the colour-octet matrix 
elements as unknown non-perturbative parameters, which is more realistic in view of 
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mcv"^ ~ A, and to exploit the constraints imposed by spin-symmetry on the leading 
contributions to these matrix elements. 
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